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Band calculations for Ce compounds with the AuCu3-type crystal structure were carried out 
on the basis of dynamical mean field theory (DMFT) . The results of applying the calculation to 
Celn3 and CeSn3 are presented as the second in a series of papers. The Kondo temperature and 
crystal-field splitting are obtained, respectively, as 190 and 390 K (CeSns), 8 and 160 K (Celn3 
under ambient pressure), and 30 and 240 K (Celn3 at a pressure of 2.75 GPa). Experimental 
results for the photoemission spectrum are reasonably well reproduced. In CeSn3, a Fermi 
surface (FS) structure similar to that obtained by a refined calculation based on the local 
density approximation (LDA) is obtained. In Celn3, the topology of the FS structure is different 
from that obtained by the LDA calculation but seems to be consistent with the results of 
de Haas-van Alphen experiments. Cyclotron mass of the correct magnitude is obtained in 
both compounds. The experimental result for the angular correlation of the electron-positron 
annihilation radiation is reasonably well reproduced on the basis of the itinerant 4/ picture. 
A band calculation for Celn3 in the antiferromagnetic state was carried out, and it was shown 
that the occupied 4/ state should have a very shallow level. 

KEYWORDS: dynamical mean field theory, band theory, de Haas van Alphen effect, ARPES, magnetic 
excitation, positron annihilation, Celn3, CeSn3 



1. Introduction 

Nonempirical band calculations for strongly correlated 
electron systems have been extensively developed on the 
basis of dynamical mean field theory (DMFT). 1 ) The 4/ 
electrons in Ce compounds are typical strongly corre- 
lated electrons. 1 ^ Recently, a DMFT band calculation 
scheme for Ce compounds was developed, which was ap- 
plied to CePd3 and CePJi3 in a previous paper cited 
as I hereafter. 2 ' In the present paper, calculations for 
Celn3 and CeSn3 will be reported. Ce compounds with 
the AuCu3-type crystal structure show a wide variety 
of 4/ electronic states from the itinerant limit to the 
localized limit. CePd3 and CeRh3 have a nonmagnetic 
Fermi liquid (FL) ground state at low temperatures. 1 ' 
Celn3 has an antiferromagnetic (AF) ground state with 
a Neel temperature of Tn = 10 K. 3 ~ 6 ) Recently, it was 
found that Tn decreases to zero at a pressure of 2.5 GPa. 
In addition, a transition to superconductivity occurs 
at Tsc = 0.2 K. 7,8 ' CeSn3 is the first heavy-fermion 
material for which an experiment on the de Haas-van 
Alphen (dHvA) effect has demonstrated the applicabil- 
ity of the itinerant 4/ band picture. 9-11 ) The Fermi sur- 
face (FS) structures of Celn3 and CeSn3 have been ex- 
tensively studied from both experimental and theoretical 
viewpoints. 12, 13 ) It will be worthwhile to check whether 
DMFT calculation gives the correct FS topology and the 
cyclotron mass obtained experimentally. 

The 4/ state splits into the j = 5/2 ground multiplet 
and the j = 7/2 excited multiplet with a separation of 
about 0.3 eV owing to the spin-orbit interaction (SOI). 
The multiplet shows crystal-field splitting (CFS) of ap- 
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proximately 100 K. The lower j = 5/2 multiplet splits 
into the (j = 5/2)r 7 doublet and the (j = 5/2)r 8 quar- 
tet in the cubic point group, which we hereafter call the 
r 7 and Tg states, respectively. It is important to take 
account of the SOI and CFS effects in 4/ systems. A 
theory named NCA/ 2 vc (noncrossing approximation in- 
cluding the f 2 state as a vertex correction) has been 
developed 14-16 - 1 and combined with the linear muffin- 
tin orbital (LMTO) method to carry out the DMFT 
band calculation. 16 ' 17 ) NCA/ 2 vc can include CFS and 
the SOI effect and also the correct exchange process of 
the f 1 — > f°, f 2 virtual excitation. The calculation gives 
an accurate order of the Kondo temperature (Tk). It was 
shown in I that the results of DMFT band calculation 
show reasonable agreement with experimental results for 
the photo emission spectrum (PES), inverse PES (IPES), 
angle-resolved PES, and inelastic magnetic excitation in 
CePd 3 and CeRh 3 . CePd 3 is a material with T K ~ 250 K, 
and CeRh-3 is one of the compounds with the most itin- 
erant 4/ states. The hybridization intensity (HI) in these 
compounds is very strong because the 4/ state hybridizes 
with a conduction band with a high partial density of 
states: the 4d states of Pd or Rh. 

In this paper, calculations for Celn3 and CeSn3 will 
be reported. The 4/ state hybridizes with bp states of 
In or Sn in these compounds. Calculated results gen- 
erally show reasonable agreement with experimental re- 
sults for the PES 18-23 ) and inelastic magnetic excitation 
by neutrons, 24 ~ 26 ) although the calculation gives some- 
what higher Tk for Celn3 than that expected from ex- 
periments, and gives a lower Tk for CeSn3 when experi- 
mental results are examined in detail. 

The FS structure at low temperatures is studied. Cal- 
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culatcd dHvA frequencies show agreement with experi- 
mental results, and cyclotron masses have the correct or- 
der of magnitude. In CeSn 3 , the calculated FS is almost 
identical to that of Hascgawa et al. 11 ^ based on the local 
density approximation (LDA). Their results show good 
agreement with those of experiments. 10 ) The FS topol- 
ogy depends on the energy levels of 4/ states relative to 
the Fermi energy (Ep). In the LDA calculation, however, 
the energy position of 4/ levels is depending on details 
of calculation. 27 ~ 32 ) For example, several LDA band cal- 
culations fail to give the topology of the experimentally 
obtained FS structure. 12 ' 13 ) The position of the 4/ lev- 
els relative to Ep is robust in DMFT because CFS is 
relatively large compared with the 4/ band dispersion. 
In Celn 3 under pressure, the calculated FS is different 
from the result of the LDA calculation 3 ' 32 ~ 37 ) A closed 
electron pocket that is centered at the L point appears, 
in contrast to the hole pocket in the LDA calculation. 3 ) 
This FS structure is, in some sense, similar to that of 
Laln 3 , 38 ~ 4 °) although this does not imply the localized 
nature of the 4/ state. The effective Kondo temperature 
markedly increases from Tk ~ 10 K at ambient pressure 
to 30 K at about p = 2.75 GPa. The topology of the 
FS at these two pressures is similar but is not identical 
even when the paramagnetic ground state is assumed at 
ambient pressure. 

The occupation number of electrons for a fixed wave 
vector k (k-ONE) is calculated for Celn 3 for comparison 
with the experimental angular correlation of the elec- 
tron positron annihilation radiation (ACAR), which has 
been considered to indicate a localized 4/ state. 41, 42 ) Ob- 
served results can also be reproduced by the 4/ band 
picture used in the present calculation. 

To examine the band structure in the AF state of 
Celn 3 , a calculation based on a LDA+U-like model is 
carried out. It is shown that a shallow level for the occu- 
pied 4/ state is needed to reproduce the dHvA frequen- 
cies. 

In §2, we briefly give the formulation on the basis of 
the LMTO method. Results of the application to Celn 3 
under pressure are shown in §3, results for Celn 3 at am- 
bient pressure are given in §4, and a calculation for the 
AF state of Celn 3 is presented in §5. The application to 
CeSn 3 is shown in §6. A summary is given in §7. 

2. Formulation 

The method of calculation is described briefly because 
its details have been given in previous papers. 2 ) We con- 
sider the excitation spectrum of the following Hamilto- 
nian: 

^ = ^ DA + |E(E^ r7 V ir7 -n^) 2 . (1) 
i r, 7 

Here, U is the Coulomb constant 43 - 1 and c^ a j p 7 is 
the annihilation operator for the atomic localized state 
dfl (r ) at site i with the 7 orbital of the T-irrcducible 

representation. The quantity n ^ A * is determined using 
the occupation number of the atomic 4/ electron per Cc 
ion in the LDA calculation. We assume that the local 



Coulomb interaction acts only on the orbital dfi , and 

the localized 4/ state </> a is approximated by the band 
center orbital </>(—). 

The excitation spectrum is expressed by introducing 
the self-energy terms, 

H D mft = Hlda+ J2 ( s r(e + W)+er-£r DA )l0I r7 >< 
i ,(r, 7 ) 

where e r is the single-electron energy level of the 4/ state 
and £p DA is the energy level in the LDA calculation. The 
self-energy Hp(e + id) is calculated by solving the auxil- 
iary impurity problem with the use of NCA/ 2 vc; its out- 
line is described in the Appendix of ref. 16. This method 
gives an accurate order of the Kondo temperature. The 
splitting of the self-energy due to the SOI and CFS ef- 
fects is considered. 

The matrix equation for the Greenian is written for a 
given wave number vector k as 

[zI-D L - DA (k)--£(z)}G(z;k) = I, (3) 

where / is the unit matrix and -Dlda (k ) is the diagonal 
matrix of the eigenenergies of Hlda with fe . The matrix 
elements of are obtained by calculating the self- 

energy operator term of eq. (2) in the manifold of the 
eigenvectors of "Hlda- 

The density of states (DOS) on the atomic 4/ state is 
given by 

p r band) (e; k ) = - ±3tr[d r G(e + iS; k )], where the pro- 
jection operator is defined as Op = Y^i 1 \ ( t ) \ , Tl \ >< 

6 s -. I. The local DOS in the DMFT band calculation 

Y i (r 7 )' 

is obtained by summing p r band ^ (e; k ) over k in the Bril- 

louin zone (BZ): p r band '(e) = J^k P^^i 6 ^ ^ )• H ere j 
N is the total number of unit cells. 

The self-consistent calculation in the DMFT is as fol- 
lows. First, we calculate the self-consistent LDA band 
by the LMTO method; potential parameters are fixed to 
those in the LDA calculation except for the / levels. (I) 
We calculate the atomic 4/ DOS ^ imp0 (e) (4/ DOS) for 
the auxiliary impurity Anderson model by the NCA/ 2 vc 
method with a trial energy dependence of the HI and 
4/ levels, then calculate the local self-energy. (II) The 
DMFT band calculation is carried out using the self- 
energy term, and the local 4/ DOS in the DMFT band is 
calculated. The calculation is iterated so that the 4/ DOS 
of the local auxiliary impurity model and the DMFT 
band satisfy the self-consistent conditions. 

The 4/ level is adjusted in the DMFT self-consistent it- 
erations under the condition that the 4/ occupation num- 
ber has a given target value, n/(rsl. target). The tempera- 
ture dependence of the Fermi energy, Ep, is neglected by 
fixing it at a value determined at a low temperature. It is 
estimated using the occupation number obtained by the 
renormalized band (RNB) calculation, in which the self- 
energy is approximated by an expansion form up to the 
linear term in the energy variable at Ep (see Appendix A 
of I for the calculation of the total occupation number). 
The target 4/ electron number n/(rsl. target) is imposed 
on the occupation number calculated directly using the 
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resolvents to stabilize the self-consistent iterations. 

3. Celn 3 under Pressure 

3. 1 Density of states 
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Fig. 1. 4/ DOS of Celn 3 with a = 8.7326 a.u. at T = 2.34 K. The 
solid line shows the total 4/ PES. The dashed line is the DOS of 
the (j = 5/2)T7 component, the dot-dash line is the DOS of the 
= 5/2)rg component, and the two-dots-dash line is the DOS 
of the j = 7/2 component. The Fermi energy E-p = 0.53231 Ry is 
indicated by the vertical dot-dash line. The inset shows spectra 
in the vicinity of Ep. 

In Fig. 1, we show the 4/ DOS at T = 2.34 K for 
Celn 3 at pressure p = 2.75 GPa with lattice constant 
a = 8.7326 a.u. 44 ) The solid line shows the total 4/ com- 
ponent of the PES (4/ PES). The dashed line is the 
DOS of the Tr component and the dot-dash line is the 
DOS of the Ts component. The two-dots-dash line is 
the DOS of the j = 7/2 component. CFS of the self- 
energy in the excited j = 7/2 multiplet is neglected. 
The vertical dot-dash line indicates the Fermi energy 
E-p — 0.53231 Ry. The inset shows spectra in the en- 
ergy region near E-p. The 4/ DOS has a large peak at 
0.557 Ry, about 0.025 Ry (0.34 eV) above E F . This peak 



Fig. 2. k-intcgratcd magnetic excitation spectrum of Celn3. Solid 
lines show the spectrum at T = 2.34 K and dashed lines show 
the spectrum in a hypothetical case where matrix elements of the 
magnetic moment are restricted within the intra-T7 manifold of 
space. Bold (thin) solid and dashed lines are calculated using the 
lattice constant a = 8.7326 a.u. (a = 8.859 a.u.). The bold (thin) 
dot-dash line is the spectrum at T = 300 K (T = 37.5 K) in the 
case of a = 8.7326 a.u. (a = 8.859 a.u.). 

mainly has the j = 7/2 character. The spin-orbit split- 
ting on the IPES side is usually enhanced in Ce systems 
as noted in I. The spectral intensity on the PES side 
mostly consists of the component. The PES has a 
sharp peak at Ep with a side peak at 0.021 Ry (0.29 eV) 
below E F . The latter corresponds to the SOI side band. 
The PES also has a broad peak at about 0.4 Ry, which 
corresponds to the atomic f 1 — > f° type excitation. The 
PES of Celn 3 at ambient pressure has been extensively 
studied in experiments. 21 ~ 23 ) A comparison with these 
experimental results will be made in a later section. 

In Fig. 2 we show the k-integrated magnetic exci- 
tation spectrum. The total magnetic excitation spec- 
trum at T = 2.34 K is shown by the bold solid line, 
which has peaks at about E = 0.0002 Ry (32 K) and 
E = 0.0015 Ry (240 K). The bold dashed line depicts 
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Tabic I. Various quantities obtained in the DMFT calculation for 
Celn 3 with a = 8.7326 a.u. at T = 2.34 K. n r lmp,) is the occupa- 
tion number in the auxiliary impurity problem of the effective HI 
of the DMFT calculation, er is the energy level in Ry. pr(^F) 



is the partial DOS at Ep in Ry 1 . Z T 



l-d?RT, T (e)/de\E F 



in the band calculation is £^ nd 
static CFS is set to be zero: e r7 



is the mass renormalization factor of the 4/ band, £p is the ef- 
fective energy of the rcnormalized band in Ry, and Tp is the 
imaginary part of the self-energy at Ep in Ry. (Tps is a very 
small positive value.) The effective energy levels are measured 
from the Fermi energy Ep = 0.53231 Ry. The spin-orbit in- 
teraction constant is £ 4 y = 7.073 X 10 — 3 Ry. The 4/ level 

0.55550 Ry. The clcctro- 
£p g . The target 4/ occu- 
pation number is n^(rsl. target) = 0.955 and the resultant oc- 
cupation number calculated using the resolvent is 0.955. The 
4/ occupation number calculated by integrating the spectrum 
is n/(intg.) = 0.967 and the obtained total band electron num- 
ber is iV (total; RNB) = 12.999. Ei no i is the characteristic en- 
ergy of the quasi-elastic excitation and BcfS is the CFS ex- 
citation energy estimated from the magnetic excitation spec- 
trum. The occupation numbers of the 4/ state in the LDA are 
n/ (Cc, LDA) = 1.061 and nj (La, LDA) = 0.107. The Coulomb 
constant U is set to be 0.51 Ry (6.9 eV). 
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the spectrum for a hypothetical case in which the ma- 
trix elements of the magnetic moment are nonzero only 
in the manifold of r 7 , and thus it may correspond to the 
excitation spectrum within the r 7 manifold. The Kondo 
temperature within the r 7 state is estimated to be 32 K 
and the CFS excitation is estimated to be 240 K. 

The parameters and the calculated values are given 
in Table I. The LMTO band parameters for states ex- 
cept for the / component arc fixed to those obtained 
by the LDA calculation. Ep is fixed to the value deter- 
mined by the occupied state in the RNB calculation. The 
occupied 4/ electron mainly has the r 7 character. The 
occupation number of the r 7 component relative to the 
T 8 component, 0.63/0.26 (2.4), is very large compared 
with 0.5 expected from the ratio of the degeneracy, but 
is smaller than the value expected from the simple model 
of CFS for an isolated ion with an excitation energy of 
240 K. The occupation of the j — 7/2 component, 0.08, 
is very small. Usually, the electrostatic potential causes 
cubic CFS in 4/ electron systems. This is neglected in 
the present calculation. The hybridization effect causes 
large CFS of 240 K. 

The target value of the occupation number on the 
atomic 4/ states, n/(rsl. target), has been tentatively set 
to be 0.955 in the present calculation. The occupation 
number on the atomic 4/ state is 1.061 in the LDA band 
calculation for Celn3. On the other hand, the occupa- 
tion number on the atomic 4/ (La) state is estimated to 
be 0.107 when we carry out the LDA calculation for a 
compound in which Ce ions are replaced by La ions. The 



Fig. 3. Hybridization intensity (HI) of Ccln3 with a = 8.7326 a.u. 
calculated by the LDA band method. The dashed line is the HI 
of the (j = 5/2)r 7 component, the dot-dash line is the HI of the 
(j = 5/2)Tg component, and the two-dots-dash line is the HI of 
the j = 7/2 component. Ep is indicated by the vertical dot-dash 
line. The inset shows HI in the vicinity of Ep. The 4/ levels are 
located at 0.71143 and 0.73620 Ry, respectively, for j = 5/2 and 
7/2. 

difference between these two 4/ occupation numbers is 
1.061 - 0.107 = 0.954. The value 0.955 is 90% of the 4/ 
occupation number (1.061) obtained by the LDA band 
calculation for Celn3 . When we perform the DMFT band 
calculation by setting n/(rsl. target) to be 0.997, 94% of 
the LDA value used for CePds, 2 ' 45 * 1 7k is estimated to 
be low (10 K). 

In Fig. 3, we show the HI obtained from the LDA 
band calculation, which is used as the starting HI in the 
DMFT calculation. It has sharp peaks. The partial den- 
sity of states on the (In)5p state has corresponding peaks, 
although they are less conspicuous. Ep is located at a 
slope of a peak as indicated in the figure. In Ccln3, the 
magnitude of the Kondo temperature is very sensitive to 
the position of the Fermi energy, although the effective 
HI in DMFT is usually smeared out and generally re- 
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duced. 46 -* We note that an instable region appears with 
9JV tot^icctron < q whcn Ep is re duccd by 0.025 Ry. There 
is a possibility that a solution with a lower Ep and an 
extremely low Kondo temperature appears, but in this 
paper we adopt a solution with higher Ep. Checking the 
relative stability of the two states and the transition be- 
tween them remains as a future work. 

3.2 RNB and wave-number-vector- dependent DOS 
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Fig. 5. Wave number vector (k ) dependence of the DOS (k-DOS) 
for Celn3 with a = 8.7326 a.u. at T = 2.34 K. k moves from the 
T point (bottom) to the X point (top) along the A axis. The 4/ 
DOS is shown. = 0.53231 Ry is indicated by the vertical dot- 
dashed line. The spectra are broadened by an imaginary factor, 
7 = 0.001 Ry, in the energy variable. The lines are plotted in 
cyclic order (solid, dashed, dot-dash, two-dots-dash, and three- 
dots-dash lines) from the Y point but the deviations from this 
regularity appear in the peaks of the spectra because some peaks 
are very large. 



Fig. 4. Band dispersions of the renormalizcd band (RNB) picture 
for Celn 3 with a = 8.7326 a.u. at T = 2.34 K. The symbols 
under the horizontal axis denote the symmetry points and axes 
of the BZ of the simple cubic (s.c.) lattice. _Bp = 0.53231 Ry is 
indicated by the horizontal dashed line. The band crossing E-p 
is the 7th band. 

In Fig. 4, we show RNB dispersions at T = 2.34 
K. The energy shift (the real part of the self-energy 
at Ep: 3?Sp(i?F)) and the mass renormalization factor 
(1 — 95RSr(e)/5e|_E F ), which are given in Table I, are 
taken into account in this calculation. The dispersions of 
the RNB can be understood as the hybridization bands 
of the 4/ and Laln 3 -like bands. 



Narrow bands with the j = 5/2 character appear 
slightly above E F . Bands with the character of j = 7/2 
appear around an energy of 0.557 Ry, which is the en- 
ergy of the j — 7/2 peak in the 4/ DOS shown in Fig. 
1, although they are not shown in Fig. 4. The 4/ bands 
located at about 0.0012 Ry above Ep mainly have the 
Tg character, and the 4/ bands located at Ep mainly 
have the T7 character. The latter is the 7th band in the 
present calculation. The state sinks below Ep near 
the T point. This is different from the result of the LDA 
calculation, in which the T7 state is located above Ep 
at the T point. 3 ' 32 ' 35 ) We note that dispersions of the 
RNB are qualitatively similar to the band structure de- 
termined by the LDA+U calculation with the shifted 4/ 
level, 47 ) but the width of the 4/ band with j = 5/2 is 
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Fig. 7. Angular dependence of the dHvA frequency of Celn3 with 
a = 8.7326 a.u. at T = 2.34K. 



Fig. 6. Fermi surface structure for the 7th band of Celn3 with 
a = 8.7326 a.u. at T = 2.34 K. The frame is -n/a < (k x , k y ) < 
w/a and —w/a < k z < 0, and the region occupied by the 7th 
band is shaded. The center and corners of the upper plane are 
the r point and M points, respectively. The center and corners 
of the lower plane are the X point and R points, respectively. 
The midpoints of the edges of the upper plane (midpoints of the 
neighboring M points) are also X points, and the midpoints of 
the edges of the lower plane (midpoints of the neighboring R 
points) are also M points. The electron pocket centered at the T 
point is closed and the electron pocket centered at the R point 
is also closed. 

about 0.005 Ry in the RNB, whereas that of the LDA+U 
band is about 0.02 Ry. 

The dispersions of the RNB correspond well to the be- 
havior of the k -dependent density of states (k-DOS) as 
already noted in I. For example, we show the 4/ compo- 
nent of the k-DOS when k moves from the T (bottom) 



Table II. Calculated and experimental dHvA frequencies and cy- 
clotron masses of Celn3 with a = 8.7326 a.u. Calculated values 
are estimated at T = 2.34 K. F is the dHvA frequency in the 
unit of 10 3 T and m c is the cyclotron mass in the unit of electron 
mass rriQ. Experimental values arc presented in parentheses, n in 
the expressions dn and a n denotes the direction of the magnetic 
field. Experimental values of oooi and dooi are estimated from 
Fig. 2 of ref. 48 at p = 2.75 GPa. Other experimental values are 
given in Table II of ref. 3 at p = 2.7 GPa. 
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to the X (top) point along the A line in Fig. 5. A peak of 
the DOS with mainly the 4/ character is located below 
Ep at the T point. Starting from the peak at the T point, 
a ridge line runs a short distance along Ep, and then sud- 
denly drops on the low-energy side because the 4/ band 
forms a hybridized band with a band of (In)5p character, 
which has strong dispersion. Near the X point, a ridge 
line with weak dispersion appears again above Ep. Note 
that we have plotted the 4/ component, not the total 
intensity. If we plot the latter, it becomes an order of 
magnitude larger than that of the T point for the wave 
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Fig. 8. Wave number vector (fe ) dependence of the occupation 
number of electrons (k-ONE) for Celn3. Solid lines are calcu- 
lated at T = 2.34 K and the bold dashed line is calculated at 
T = 60 K in the case of a = 8.7326 a.u. The bold solid line is 
the total occupation number and the thin solid line is the oc- 
cupation number of the 4/ component. The thin dashed line is 
calculated at T = 60 K for Celn3 with a = 8.859 a.u. For the to- 
tal occupation number, 13 is subtracted and the scale on the left 
side is used. For the 4/ occupation number, 1 is subtracted and 
the scale on the right side is used. In the calculation of the total 
k-ONE, an extra broadening factor of 7 = 0.0075 Ry is added in 
the k-DOS to avoid the numerical error in the integration of the 
sharp spectrum. The k-ONE of the 4/ component has a small 
random oscillation due to the numerical noise because the extra 
broadening is not added. 

number vector k ~ 0.6 x ^ (the wave number vector 
corresponding to the 3rd solid line from the bottom) in 
the vicinity of Ep because the contribution of the (In)5p 
component is large. The spectral intensity is greatly re- 
duced to about Z ~ 1/30 when the 4/ component is 
dominant. 

In Fig. 6, we show the FS structure of the 7th band cal- 
culated from the RNB dispersion of DMFT at T=2.34K. 
Two electron pockets appear, a smaller one centered at 
the T point and larger one centered at the R point. A 



cagelike connected-electron sheet exists in the space be- 
tween the two electron pockets. A similar FS structure 
has been obtained from the LDA+U calculation with the 
shifted 4/ level. 47 ) 

The small electron pocket at the T point is closed, 
i.e., it does not touch other sheets. On the other hand, 
the pocket touches the cagelike sheet in the LDA. In 
0Q addition, it has a hollow (a hole pocket) inside it in the 
LDA calculation because the Tr state is located above Ep 
at the r point. 3 ) The large electron pocket at the R point 
is closed in both the DMFT and LDA calculations. 3 ' 

In Fig. 7, we show the dHvA frequency calculated 
' 50 from the RNB dispersion in the DMFT calculation. The 
electron pocket at the T point gives signals at about 
4 x 10 3 Tcsla (T) and the large electron pocket cen- 
tered at the R point gives signals at about 1 x 10 4 T. 
0Q These signals correspond to the d and a branches in ex- 
perimental results. 3 ' 48 ' The dHvA frequency (F) and cy- 
clotron mass (m c ) for magnetic fields along the princi- 
pal axes are listed in Table II. The calculated frequen- 
cies show good agreement with experimental results, al- 
osfihough they are somewhat larger. The agreement of the 
cyclotron masses is less good, but the calculation gives 
a comparable magnitude. Many other branches appear 
as seen in Fig. 7. For example, in the case of H || (001) 
we have a branch with F = 1.4 x 10 3 (m c = 63mo), 
two branches with F = 0.7 x 10 3 (m c = 12m , 16m ), a 
branch with F = 0.5 x 10 3 (m c = 7 mo), and a branch 
with F — 0.4 x 10 3 (m c = 16m ). The curvature fac- 
tor (the second derivative of the cross-section area) of 
these branches is not so large. However, the calculated 
results are less definite because they depend strongly on 
fine structures of the cagelike FS. 

Recently, the momentum density of occupied electrons 
has been extensively studied by the ACAR. Biasini and 
co-workers have concluded that the momentum density 
of Celn3 in the paramagnetic state coincides with the 
predictions of band calculations in which the 4/ elec- 
trons are regarded as localized. 41,42 ) For example, let us 
imagine unhybridized Laln3-like bands by hypothetically 
removing the 4/ components from Fig. 4. When k moves 
from T to X along the A line, a conduction band cuts 
the Fermi energy from downward to upward, then an- 
other band cuts it from upward to downward. A similar 
band dispersion appears around the T point. Therefore, a 
small electron pocket appears with a characteristic shell- 
like hole region around it in the Laln3-like band. 3 ' 12 ' 49 ) 
A valley of momentum density appears halfway between 
r and X points. Such a valley is observed in the experi- 
mental results of ACAR. 42 ) On the other hand the LDA 
band of Celn3 shows a sharp increase halfway from the 
r point to the X point, or to the R point because the 
electron pocket at the T point has the hollow of the hole 
pocket inside it. 3 ) 

In Fig. 8, we show the total k-ONE on the symmetry 
points and axes obtained by the DMFT band calculation 
by a bold solid line. We note that k-ONE has a valley 
when k moves from the T point to the X point. In the 
dispersions shown in Fig. 4, the T point is included in the 
electron pocket and the X point is included in the hole 
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region. Therefore, k-ONE would sharply decrease at the 
point where the band crosses the Fermi energy when k 
moves from the T point to the X point along the A line 
if the band dispersions in Fig. 4 were dispersions of the 
noninteracting bands. The k-ONE given in Fig. 8 does 
not always show variations corresponding to the inter- 
section of the RNB band with the Fermi energy. This is 
partly ascribed to the large mass renormalization factor 
of 4/ bands: the discontinuity of k-ONE is reduced to 
Z ~ 1/30 when the crossing hybridized bands mainly 
have the 4/ character. k-ONE shows a steep but nonsin- 
gular variation when a Laln 3 -like band cuts the Fermi 
energy. Another factor may be the extra broadening in- 
troduced to avoid a numerical error in the integration 
of the sharp spectrum of the k-DOS. In Fig. 8, we also 
show the variation in k-ONE of the partial k-DOS of 
the 4/ component by the thin solid line. In this calcu- 
lation, only the lifetime broadening of the DMFT band 
state is considered, and the extra broadening factor is 
not added because the 4/ k-DOS has a larger width. We 
can recognize a relatively small but rather rapid vari- 
ation of k-ONE corresponding to the Fermi function. 50 ) 
We may conclude that the main variation of k-ONE is de- 
termined by the dispersion of the hypothetical Laln 3 -like 
band. The calculated result in DMFT does not appear to 
contradict the experimental results of ACAR. We have 
also calculated k-ONE at higher temperatures and also 
at ambient pressures. The total k-ONE in these cases is 
similar to the result shown by the bold solid line in Fig. 8, 
except that the rapid variation in the 4/ component is 
smeared. 

We have noted that the calculated dHvA frequency 
is consistent with experimental results. The dHvA fre- 
quency reflects the singular change in the momentum 
density, thus it can detect the FS even when the mass 
renormalization factor is very large. An ACAR experi- 
ment with a very fine resolution at very low tempera- 
tures should be carried out. The calculation of ACAR 
including actual momentum density remains as a future 
work. 

When we calculate the RNB dispersion at T = 150 
K, the 4/ band shifts slightly to the high-energy side. 
This shift of the 4/ band with increasing temperature 
has already been noted in I. The 4/ state at the V point 
shifts to above Ep. Therefore, a hole pocket appears in- 
side the electron pocket centered at the T point. The 
primary structure of the FS of the RNB at T = 150 K is 
similar to that of the LDA band of Celn3, not to that of 
Laln3, except for some differences in the fine topology. 
Note that the RNB picture has only limited meaning at 
high temperatures because the imaginary part is large. 
However, the trace of broad peaks in the k-DOS shows a 
shift corresponding to that of the RNB dispersion as al- 
ready shown for CePd3 in I. The calculated k-ONE does 
not show any indication of the hole pocket at the T point 
even at high temperatures, as already noted. 

4. Celri3 at Ambient Pressure 

In this section we show calculated results for Celn3 
with the lattice constant a = 8.859 a.u. under ambient 
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Fig. 9. 4/ DOS of Celn 3 with a = 8.859 a.u. at T = 100 K. For 
the definition of lines, see the caption of Fig. 1. The Fermi energy 
-Ef = 0.50088 Ry is indicated by the vertical dot-dash line. 



pressure. 51 ^ Celn 3 enters the AF state at Tn = 10 K at 
ambient pressure, 4 ) but we tentatively assume the para- 
magnetic state in the present section. In Fig. 9, we show 
the 4/ DOS calculated at T = 100 K. It shows peaks 
of the Kondo resonance and the SOI side band near the 
Fermi edge and a broad peak at about 0.34 Ry, whose 
binding energy is about 0.17 Ry (2.3 eV) from Ef. Kim 
and co-workers have observed the PES at the 4c?-4/ and 
3<i-4/ resonance thresholds and carefully separated the 
bulk and surface components. 21, 22 ) The bulk component 
of the PES has a broad peak with a binding energy 2 eV 
from Ef- The calculated intensity of the Kondo reso- 
nance part relative to the SOI side band peak is larger 
than the experimental value. The magnetic excitation 
spectrum is shown in Fig. 2 by thin lines. CFS is esti- 
mated to be 160K (~ 0.001 Ry - 14 meV) from the peak 
position of the thin solid line, and Xk is estimated to be 
8 K (~ 0.00005 Ry - 0.7 meV) from the peak position of 
the thin dashed line. In an experiment on neutron scat- 
tering at T = 20 K, a broad peak of CFS excitation is 
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Tabic III. Various quantities obtained in the DMFT calculation 
for Ccln 3 with a = 8.859 a.u. at T = 2.34 K. For the defini- 
tion of notation, see the caption of Table I. The Fermi energy 
is E F = 0.50088 Ry and AEg = K. The spin-orbit interac- 
tion constant is ( 4 / = 7.057 X 10~ 3 Ry. The 4/ level in the 
band calculation is e\^ nd = 0.52212 Ry. The target 4/ occu- 
pation number is n/(rsl. target) = 0.969 and the resultant oc- 
cupation number calculated using the resolvent is 0.969. The 
4/ electron number calculated by integrating the spectrum is 
ri/(intg.) = 0.977 and the calculated total band electron num- 
ber is N (total; RNB) = 13.000. 71/ (Ce, LDA) = 1.076 and 
n/ (La, LDA) = 0.104. The Coulomb constant U is set to be 
0.51 Ry (6.9 eV). Pps is a very small positive value. 
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observed at about 13 meV, and a quasi-elastic excitation 
peak with a width of about 2 meV is also observed. 26 ) 
In the calculation at T = 37.5 K shown by the thin dot- 
dash line, the width of the quasi-elastic excitation peak 
is about 0.0002 Ry (2.7 mcV). The calculated value has 
the correct magnitude, but Tk seems to be somewhat 
higher than the experimental value. 

We have carried out a DMFT band calculation at 
T = 2.34 K assuming the paramagnetic state, although 
Celn3 enters the AF state at low temperatures. The pa- 
rameters and calculated values are given in Table III. 
The mass enhancement factor is small compared with 
that when p = 2.75 GPa, although the ratio of Tk is 
large, 32 K/8 K - 4. 

The overall features of the RNB band dispersion are 
almost identical to those when p = 2.75 GPa, which are 
shown in Fig. 4, but the band is located relatively 
deep from Ep when p = 0. The state is located below 
Ep at the T point, and a closed small electron sheet cen- 
tered at the r point also appears under ambient pressure. 
One important difference is that the band is located 
below Ep at the M point. Therefore, the large electron 
sheets centered at the R-point touch each other at the 
necks at M points. The FS structure has been presented 
in Fig. 5 of ref. 52. dHvA signals with frequencies of ap- 
proximately 3.8 x 10 3 T and cyclotron masses m c ~ 15too 
appear, which correspond to orbits on the small electron 
pocket centered at the T point. This will be discussed in 
the next section. 

In Fig. 8, we show k-ONE at T = 60 K under am- 
bient pressure by the thin dashed line. It has a valley 
when k moves from the T point to the X point. When k 
moves from the T point to the R point, it shows a shal- 
low minimum and a steep increase. The calculated result 
for total momentum density in DMFT does not contra- 
dict the experimental results of ACAR 41,42 ) without the 
exceptional assumption of the localized 4/ state. 
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Fig. 10. Band dispersions for the AF state of Celn3 with a = 
8.859 a.u. Ep = 0.50610 Ry is indicated by the horizontal dashed 
line. Dispersions arc depicted using the symmetry points and 
axes of the BZ of the s.c. lattice to illustrate the zone-folding 
effect in the AF state. The bands crossing Ep are the 13th and 
14th bands. 



5. Celn 3 in AF State 

At ambient pressure, Celn3 undergoes a transition to 
the AF state at Tn = 10 K. In this section, the band cal- 
culation assuming the AF state is carried out. The order- 
ing of Ti orbitals with Q = (1/2, 1/2, 1/2) is assumed. 
The orbital of the state with the polarization paral- 
lel to the (1,1,1) direction has low energy in a sublattice, 
and in another sublattice the orbital with the inverse po- 
larization has low energy. 4 ^ The energy level of these oc- 
cupied orbitals is chosen to be E^ ^ — Ep — —0.0042 Ry, 
where Ep = 0.5061 Ry. The energy level of unoccupied 
states is set to be a higher energy of 0.2058 Ry above Ep. 
To see the characteristic features of the band structure 
in the AF state, a calculation similar to the LDA+U 
calculation is carried out using these energy levels but 
without the self-consistent condition on U-tcrms. 53 ) We 
note that if the energy level of whole 4/ orbitals is set 
to be that of the unoccupied state, dispersions of the 
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Fig. 11. Fermi surface structure for the 13th band of the AF state 
of Celn3 with a = 8.859 a.u. The frame is —n/a < (k x , ky) < n/a 
and —n/a < k z < 0, and the hole region is shaded. The center of 
the upper plane is the T point and the corners of the lower plane 
are R points. These two points are equivalent in the AF state. 
For the names of axes and points, see the caption of Fig. 6. 

Laln 3 -like band are obtained. The result in the AF state 
is shown in Fig. 10. It is depicted in the frame of the 
BZ of the original simple s.c. lattice, therefore, the fold- 
ing of bands due to the AF ordering can be seen. 54 ) For 
example, the center of the A axis is located on the zone 
boundary of the BZ in the AF state, and the R point and 
T point are equivalent. Band dispersions are symmetric 
with respect to the central point on the A axis. 55 -* 

Note that we have chosen a very shallow energy level 
for the occupied 4/ states. When we carry out the cal- 
culation assuming deeper levels for occupied states, rec- 



Fig. 12. Fermi surface structure for the 14th band of the AF state 
of Celn3 with a = 8.859 a.u. The frame is —n/a < (k x , k y ) < n/a 
and —n/a < k z < 0, and the occupied region is shaded. For the 
names of axes and points, see the caption of Fig. 11. 

ognizable changes in band structures from those of the 
Laln 3 -like band do not appear in the energy region near 
Ep, although new zone boundaries of the AF state are de- 
fined in a mathematical sense. The hybridization effect in 
the vicinity of Ep is very small. When we choose shallow 
4/ levels for the occupied Tr orbitals, actual changes in 
band structures appear in the energy region near Ep . The 
bands showing a steep increase above Ep at the midpoint 
of the A axis originate from the folding of a band, which 
gives the large electron pocket centered at the R point 
in the Laln 3 -like band. These steep bands with non-4/ 
character and the flat band with the character slightly 
below Ep show typical band repulsion behavior due to 



J. Phys. Soc. Jpn. 



Full Paper 



Osamu Sakai, Hisatomo Harima 



11 



10 



4 CeIn 3 AF a=8.859au 



o 

CD 



** 1 




■*«*»+ 




[101] [001] [111] [110] 

Field Angle 



Fig. 13. Angular dependence of the dHvA frequencies for the AF 
state of Celn3 with a = 8.859 a.u. 



the c — f hybridization. On the other hand, the flat band 
does not show appreciable band repulsion with the band 
that forms the surface of the small electron pocket cen- 
tered at the T point, as seen on the S and A axes in 
Fig. 10. This small electron pocket corresponds to the 
small electron pocket at the T point in the Laln3-likc 
band. The latter has a shell-like hole sheet around it. 3 ) 
The band forming the outer surface of the shell-like hole 
shows band repulsion due to the c — / hybridization. The 
hole region extends in the outward direction and reaches 
the X point. The band dispersion of the occupied 4/ state 
is very similar to the band in Fig. 4, especially for the 
region forming the small electron pocket centered at the 
T point. 

In Figs. 11 and 12, FSs are presented. Let us first con- 
sider the folding of the Laln3-like band in the BZ of the 
AF state. The small electron pocket centered at the T 
point with the shell-like hole sheet around it (see Fig. 
16(a) of ref. 3) is mapped on the R point. This hole 
sheet has a window (electron region) in the (111) direc- 
tion. The large electron pocket centered at the R point 
(see Fig. 16(b) of ref. 3) is also mapped on the T point. 
The large pocket includes almost the entire shell-like hole 
sheet inside it. Next we consider the c — / hybridization 
effect by introducing the shallow band. The surface 
of the small electron pocket is hardly changed but the 
outer surface of the shell-like hole sheet moves in the 
outward direction away from the small electron pocket, 



i.e., the hole region becomes wide. In the (111) direction, 
the window of the shell-like hole changes to a hump of 
the small electron pocket. 56 ) The small electron pocket 
is the FS of the 14th band shown in Fig. 12. 

The surface of the large electron pocket disappears on 
the A axis because the corresponding 14th band is shifted 
up above Ep as seen in Fig. 10. The outer surface of the 
shell-like hole sheet also disappears on the A axis since 
the 13th band is located below Ep. The 13th band rises 
up near Ep, but is still located slightly below Ep along 
the X and S axes. However, unoccupied regions usually 
exist in the 13th band in general directions. Such a hole 
region can be seen on the Z axis in Fig. 10. The hole sheet 
of the 13th band is shown in Fig. 12. It has large "holes" 
(electron regions) along the (111) direction because the 
entire 13th band is located below Ep on the A axis. It 
also has "holes" on the A and S axes. 

In Fig. 13, we show the angle dependence of the dHvA 
signals, which is calculated neglecting the rotation of the 
sublattice polarization under the applied magnetic field. 
The branch with F = 4 x 10 3 T comprises the signals 
due to the small electron pocket (14th band), and it 
corresponds to the d branch in experiments. 3 ' 12 . 51 . 57 - 61 ) 
The experimental dHvA frequencies (cyclotron masses) 
of the d branch are 3.2 x 10 3 T (2.9m ) for the field di- 
rection H || (001), 3.6 x 10 3 T (16m ) for H || (110), 
and 2.9 x 10 3 T (3.0m ) and 3.1 x 10 3 T (13m ) {d? 
branch) for H \\ (111). 3 ' The calculated values are, re- 
spectively, 3.6 x 10 3 T (0.54to ), 4.1 x 10 3 T (1.4m ), and 
3.2 x 10 3 T (0.52to ) and 3.5 x 10 3 T (1.2m ). The calcu- 
lated frequencies show a reasonably good correspondence 
to the experimental values, although they are somewhat 
larger. The calculated masses are smaller than the ex- 
perimental values because the renormalization due to the 
correlation effects is not included. We note that the small 
electron pocket centered at the T point in the paramag- 
netic state also gives dHvA signals similar to that of the 
d branch in the AF state. The calculated values are, re- 
spectively, 3.7 x 10 3 T (11.5m ), 4.0 x 10 3 T (37.7to ), 
and 3.3 x 10 3 T (12.1m ) and 3.5 x 10 3 T (28.5m ). The 
cyclotron masses estimated in the paramagnetic state are 
3 — 5 times larger than experimental values. For the field 
direction where the experimental values of the mass are 
large, the calculated masses are also relatively large in 
both the AF and paramagnetic states. 

The branch with F ~ 1.5 x 10 2 T for H \\ (001) in 
Fig. 13 has similar angle dependence to the s branch 
named by Endo et al. 59 ^ This branch runs on arms that 
connect sheets with a cross-section area having the ap- 
pearance of an upside-down spade on the k x = ir/a 
plane in Fig. 11 if the field is applied along the x di- 
rection. The centers of the orbits for H || (001) are lo- 
cated at (0.10,0.12,0.38) and (0.10,0.38,0.12) and also at 
equivalent points in the BZ. Many branches appear with 
dHvA frequencies of approximately 8 x 10 2 and 5 x 10 T. 
Branches with similar dHvA frequencies have been ob- 
served in experiments. 12 ' 59 ~ 61 ) However, their calculated 
angle dependence of them does not appear to match ob- 
served ones. In the present calculation, we have not found 
signs of hole pockets with the shape of oblate ellipsoids, 
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which were predicted in ref. 61 

The renormalization effect has been not included in 
the calculation in the present section. A very shallow 4/ 
level is needed to realize the AF effect on the band struc- 
ture. This indicates that the DMFT calculation must 
even be performed for the AF state in a future work. 
The FS structure corresponding to the d branch should 
not change drastically, even in such a calculation, be- 
cause the c — / hybridization has a weak effect on it. On 
the other hand, the FS structure of the 14th band will 
change. 

6. CeSn 3 
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Fig. 15. k-integrated magnetic excitation spectrum of CeSn3. 
Solid line shows the spectrum at T = 18.75 K and the dot-dash 
lines shows the spectrum in a hypothetical case where matrix ele- 
ments of the magnetic moment are restricted within the intra-T7 
manifold of space. The dashed line is the spectrum at T = 300 K. 



Fig. 14. 4/ DOS of CcSn 3 at T = 18.75 K. For the definition of 
lines, see the caption of Fig. 1. = 0.53160 Ry is indicated by 
the vertical dot-dash line. 

In Figs. 14 and 15, we show the 4/ DOS and mag- 
netic excitation spectra calculated for CeSn 3 at T = 
18.75 K, respectively. In the magnetic excitation spec- 
trum, the quasi-elastic component due to the Kondo ef- 
fect and the CFS component merge into a broad sin- 
gle peak at about 0.0025 Ry (34 meV). However, a 
weak swelling due to the quasi-elastic excitation in the 
Ti levels can be seen. We estimate Tk and CFS, re- 
spectively, to be 190 K (16 meV - 0.0012 Ry) and 
390 K (34 meV - 0.0025 Ry) from the peak positions 
of the dot-dash line and solid line, respectively. A broad 
peak was observed at about 40 meV in an experiment on 
inelastic neutron scattering at T=20 K. 24 ~ 26 ) A shoul- 



der like structure was also observed at approximately 
15mcV. The calculated values show good correspondence 
with the experimental results. The calculated value of 
Tk is very sensitive to the lattice constant. If we use 
a = 8.4515 a.u, which was used in the calculation of 
Hasegawa et al, Tk is estimated to be higher than 10 3 K. 
Instead, we use a = 8.921 a.u. which was determined by 
Umehara et al., 10 ^ and obtained the correct magnitude 
of T K . 

The 4/ DOS has a sharp peak at Ep and a small 
peak due to the SOI side band. A peak at a deep en- 
ergy of 0.27 Ry, whose binding energy from Ep is about 
0.26 Ry (3.5 eV) also appears. At first glance, the cal- 
culated spectrum is very similar to the spectrum under 
the 4d — 4/ resonance condition, but the experimental 
spectrum is considered to be a superposition of the sur- 
face and bulk contributions. 20 ~ 23 ) The calculated spec- 
tral shape around Ep has the characteristics of the bulk 
component. On the other hand, the binding energy of 
the deep peak is large compared with the experimentally 
observed value of about 2.5eV. 21 ^ This is mainly due to 
the fact that we chose the target value of the 4/ occupa- 
tion number to be n/(rsl. target) = 0.98. The 4/ occupa- 
tion numbers in CeSn 3 and LaSn 3 are, respectively, 1.184 
and 0.167 in the LDA calculation. The difference between 
them is 1.017. If we use n/(rsl. target) = 1.02, the bind- 
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Fig. 16. Band dispersions in the renormalizcd band (RNB) pic- 
ture for CcSn 3 at T = 18.75 K. E F = 0.53160 Ry is indicated by 
the horizontal dashed line. The 4-fold (2-fold) degenerate state 
at about 0.5334 Ry (0.5312 Ry) at the F point is the T s (r 7 ) 
state. The bands crossing Ep are the 8th and 9th bands. 

ing energy from Ep becomes 4.2 eV because we need a 
deep 4/ level to increase the 4/ occupation beyond l.O. 45 * 1 
In this case, T K and CFS are respectively estimated to 
be 110 and 280 K from the magnetic excitation spec- 
trum. When we choose n/(rsl. target) = 0.96, the bind- 
ing energy of the deep peak becomes 0.23 Ry (3.1 eV). 
In this case, Ik and CFS are respectively estimated to 
be 280 and 480 K. However, the shoulder structure due 
to the SOI side band 20 ' 21 - 1 becomes invisible. We have 
tentatively used n/(rsl. target) = 0.98 in this paper. The 
intensity ratio of the f 1 peak to the total / spectral in- 
tensity in the IPES is obtained to be 0.12, whereas a 
value of 0.18 was obtained experimentally. 19 ' The pa- 
rameters and calculated values in the DMFT calculation 
for CeSn 3 are given in Table IV. The occupation of Y$ is 
greater than that in the case of Celn 3 , although the CFS 
excitation energy is large. 

In Fig. 16, the RNB is shown. Band dispersions near 
the Fermi energy are similar to those of Hasegawa et 



Table IV. Various quantities obtained in the DMFT calculation 
for CcSn 3 with a = 8.921 a.u. at T = 18.75 K. For the def- 
inition of notation, see the caption of Table I. The Fermi en- 
ergy is E F = 0.53160 Ry and AE$ = K. The spin-orbit 
interaction constant is £4/ = 7.008 X 10~ 3 Ry. The 4/ level 
in the band calculation is e^ nd = 0.54968 Ry. The target 
4/ occupation number is nj (rsl. target) = 0.98 and the resul- 
tant occupation number calculated using the resolvent is 0.980. 
The 4/ electron number calculated by integrating the spectrum 
is nj(intg.) = 0.990 and the calculated total band electron 
number is iV (total; RNB) = 15.995. n f (Ce, LDA) = 1.184 and 
(La, LDA) = 0.167. The Coulomb constant U is set to be 
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0.382 


0.142 


er(Ry) 


-0.15825 


-0.16101 


-0.18106 


Pr (E F ) (Ry- 1 ) 


23.9 


10.6 


1.1 




14.7 


12.5 


3.3 


er(Ry) 


0.0204 


0.0422 


0.1763 


fr(Ry) 


0.78 x 10-' 


1 0.0 x 10~ 3 


0.16 x 10" 3 


£ inc l = 16 meV, 


Ecfs = 


34 meV 





Table V. Calculated and experimental dHvA frequencies and cy- 
clotron masses of CeSn3. Calculated values are estimated at 
T = 18.75 K. F is the dHvA frequency in the unit of 10 3 T 
and m c is the cyclotron mass in the unit of electron mass mo. 
Experimental values are presented in parentheses. The labeling 
of branches follows that of ref. 11. The experimental values of 
dHvA frequency are taken from ref. 10 and the cyclotron masses 
arc taken from ref. 12. 



branch F m c 



aooi 


8.4 


13.6 


1110 


8.4 (8.9) 


9.3 (4.2) 


am 


7.8 (8.0) 


7.5 (3.8) 


bmn 


10.3 (9.9) 


4.1 (2.7) 


C110 


7.4 (7.4) 


12.7 (6.3) 


dm 


6.4 (6.2) 


11.2 (6.3) 



al. based on the LDA calculation. n ) The Y$ and Ty 
states are located above E-p at the Y point in their re- 
sult. The fact that the Tg state is located above Ep leads 
to the characteristic features of an electron sheet of the 
9th band of CeSn 3 (see Fig. 5 of ref. 11). The sheet is 
centered at the Y point and is essentially spherical but 
deeply concave in the (111) directions. The eight con- 
cavities are connected with each other at the Y point 
through a small hollow when the Yg state appears above 
Ep. This result is robust in the present DMFT calcula- 
tion because CFS is relatively large compared with the 
dispersion of the 4/ band. On the other hand, a very 
careful calculation is needed to obtain this result in the 
LDA. It is known that the result of Hasegawa et al. is con- 
sistent with experimental results. 10, 12, 13 - ) In the present 
calculation, the very small hole pocket (the p branch in 
ref. 11) that originates from the 7th band and is centered 
at the r point does not appear because the Y? state is 
located below Ep, in contrast to the result of the LDA 
calculation. 

The dHvA frequencies and cyclotron masses are given 
in Table V. Branch a originates from the 8th band which 
gives a hole pocket at the R point. The calculated dHvA 
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frequencies show agreement with experimental values, 
but the cyclotron masses are about twice the experimen- 
tal values. The 8th band gives the dHvA signal of a oi in 
the calculation, but it has not been observed in experi- 
ments. The reason for this disappearance may be partly 
ascribed to the large cyclotron mass and large curvature 
factor as noted in ref. 11. However, the mass is not par- 
ticularly large compared with that of the Cuo signal, and 
the curvature factor of about 1.3 times of cuo is also not 
particularly large. At present, the reason for this disap- 
pearance is not clear. 

Many branches appear that originate from the 9th 
band. Branches with larger dHvA frequency show cor- 
respondence with experimentally observed branches as 
listed in Table V, but it becomes difficult to find experi- 
mental candidates for low-frequency ones because the de- 
viation of the frequency and angle dependence becomes 
conspicuous. 

7. Summary and Discussion 

We have studied the electronic structures of Ccln3 and 
CeSn3 on the basis of the DMFT calculation. The Kondo 
temperatures and CFS energies arc estimated to be, re- 
spectively, 8 and 160 K for Celn 3 at ambient pressure, 
and 190 and 390 K for CeSn3 from the k-integrated mag- 
netic excitation spectrum. These values agree approxi- 
mately with experimental values, but they are somewhat 
higher for Celn3 and somewhat lower for CeSn3. The 
magnetic excitation spectrum of CeSn3 shows a single 
broad peak, but this is a merged peak consisting of the 
CFS component and the quasi-elastic component due to 
the Kondo effect. At a pressure of p = 2.75 GPa, for 
which the volume contraction AV/V is about —0.04, the 
Tk and CFS of Celri3 markedly increase to 32 and 240 K, 
respectively. The decrease in the mass enhancement fac- 
tor to about 28/42 ~ 0.7 is, however, not particularly 
large. It is noted that the Kondo temperature of Ccln3 
is very sensitive to the position of Ep, because the HI 
has a very sharp peak, and Ep is located at a slope of 
a peak. There is a possibility that a transition between 
states with different Tk appears. 

In the PES of Celn 3 at ambient pressure, a broad peak 
corresponding to the /° final state appears at a binding 
energy 2.3 eV from E-p, consistent with the experimen- 
tal value of about 2 eV. The intensity of the Kondo 
resonance part relative to the SOI side band peak is 
somewhat larger than that obtained experimentally. The 
calculated Kondo temperature is somewhat higher than 
that expected from the PES experiment. 

The calculated binding energy of the f° peak of CeSn 3 
is 3.5 eV from E-p and is large compared with the experi- 
mental value of 2.5 eV . When we carry out a calculation 
assuming a shallower 4/ level to reproduce the experi- 
mental binding energy, the shoulder structure of the SOI 
side band disappears. The HI in the present calculation 
is larger than that expected from experiments. 

In Celn3, the DMFT band calculation gives FS struc- 
tures different from those obtained by the LDA cal- 
culation, whereas in CeSn3 it gives an almost equiva- 
lent FS to that in the LDA calculation. At the pres- 



sure p = 2.75 GPa, Celri3 has two closed electron pock- 
ets, which are centered at the R point and T point. 
These respectively lead to the a branch and d branch 
observed in experiments. Their calculated dHvA frequen- 
cies show correspondence with experimental results, and 
the cyclotron masses also have the correct magnitude. In 
the spaces between these two electron pockets, a cage- 
like electron sheet appears with a complicated connected 
structure. Many branches corresponding to the orbits on 
this FS are expected, although they have not yet been 
observed in experiments. The energy level of the 4/ state 
gradually increases relative to Ep with increasing pres- 
sure, and thus the Ty state will rise up above Ep at very 
high pressures. In this case, the topology of the FS will 
become similar to that in the LDA calculation. 

The occupation number of electrons for a fixed wave 
vector k (k-ONE) of Celn 3 was calculated for compari- 
son with the experimental result of ACAR. When a band 
that mainly has the 4/ character cuts the Fermi energy, 
the sharp change in k-ONE is smeared out at tempera- 
tures T of about 60 K. This is because the mass renor- 
malization factor is large and the lifetime broadening is 
also large. In such cases, k-ONE mainly reflects that of 
Laln 3 -like bands. The calculated results do not appear 
to be not contradict the experimental results. The exper- 
imental result of ACAR does not necessarily imply the 
localization of the 4/ state in Celn 3 at ambient pressure. 

The characteristic shape of the FS of the 9th band in 
CeSn 3 is reproduced by the present DMFT calculation 
because the r 8 state is located above Ep at the T point. 
This type of band dispersion naturally occurs in DMFT 
since CFS is relatively large compared with the 4/ band 
width. The calculated dHvA frequencies show agreement 
with the results of experiments for branches with larger 
frequency. The calculated cyclotron masses are about 
twice the experimental values for these branches. It is 
difficult to find a one-to-one correspondence of signals 
for lower-frequency branches. 

Band dispersions for the AF state of Celn 3 were exam- 
ined on the basis of a LDA+U-like calculation. A shal- 
low energy level for the occupied state is needed to cause 
a sizable change of FS structures of the AF state and 
those of the paramagnetic state. A small electron pocket 
centered at the T point (and also at the R point) gives 
dHvA signals corresponding to the d branch observed in 
experiments. Calculated dHvA frequencies show agree- 
ment with those obtained experimentally, but the cal- 
culated cyclotron masses are small. This pocket origi- 
nates from the small electron pocket with a shell-like 
hole sheet around it in the Laln3-like band structure. 
The small pocket centered at the T point in the DMFT 
band of the paramagnetic state gives almost equal dHvA 
frequencies to those of the AF state. The calculated cy- 
clotron masses in the paramagnetic state are about 3 — 5 
times larger than the experimental values. The masses 
are enhanced even in the AF state, but the enhancement 
factor is considerably reduced compared with that in the 
paramagnetic state. 

Let us summarize the results of the DMFT band cal- 
culation for AuCu3-type Ce compounds based on the 
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LMTO+NCAfVc method. The general features of the 
experimental results are reproduced by the calculations. 
The FS structures are almost identical to those obtained 
by the LDA calculation for CePd3, CeRh 3 , and CeSn 3 . 
On the other hand, the FS structure is different from that 
in Celn 3 . The result seems to be consistent with those of 
experiments. Experimental findings that have been con- 
sidered previously to indicate the localized 4/ state are 
reproduced within the framework of the 4/ band picture. 
The present DMFT calculation gives a higher Kondo 
temperature than that obtained from the detailed anal- 
ysis of experiments for CePd 3 and CeRh 3 , 62 ) a lower one 
for CeSn 3 , and a considerably higher one for Celn 3 . 

There is arbitrariness in the choice of the target value 
of the 4/ occupation number, n/(rsl. target). This in- 
duces an uncertainty in the calculated results. Even when 
we examine calculations by changing this value within a 
range, it is not easy to reproduce various experimental 
results consistently and quantitatively. A major origin of 
the discrepancy may be that the calculated HI is slightly 
stronger than that expected from the experimental re- 
sult. Paying attention to these features, further applica- 
tions to systems with more complex crystal structures 
are desirable. Calculations of the 4/ band state of Cc 
compounds in the AF state will be carried out in the 
near future. 
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